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ABSTRACT 

The detection of ionized bubbles around quasars in redshifted 21 -cm maps is possibly one 
of the most direct future probes of reionization. We consider two models for the growth of 
spherical ionized bubbles to study the apparent shapes of the bubbles in redshifted 21 -cm 
maps, taking into account the finite light travel time (FLTT) across the bubble. In both models 
the bubble has a period of rapid growth beyond which its radius either saturates or grows 
slowly. We find that the FLTT, whose effect is particularly pronounced for large bubbles, 
causes the bubble's image to continue to grow well after its actual growth is over There are 
two distinct FLTT distortions in the bubble's image: (i) its apparent center is shifted along the 
line of sight (LOS) towards the observer from the quasar; (ii) it is anisotropic along the LOS. 
The bubble initially appears elongated along the LOS. This is reversed in the later stages of 
growth where the bubble appears compressed. 

The FLTT distortions are expected to have an impact on matched filter bubble detection 
where it is most convenient to use a spherical template for the filter We find that the best 
matched spherical filter gives a reasonably good estimate of the size and the shift in the center 
of the anisotropic image. The mismatch between the spherical filter and the anisotropic image 
causes a degradation in the SNR relative to that of a spherical bubble. The degradation is in 
the range 10 — 20% during the period of rapid growth when the image appears elongated, and 
is less than 10% in the later stages when the image appears compressed. 

We conclude that a spherical filter is adequate for bubble detection. The FLTT distortions 
do not affect the lower limits for bubble detection with 1000 hr of GMRT observations. The 
smallest spherical filter for which a detection is possible has comoving radii 24Mpc and 
33 Mpc for a 3(7 and 5cr detection respectively, assuming a neutral fraction 0.6 at z ~ 8. 

Key words: methods: data analysis - cosmology: theory: - diffuse radiation 



1 INTRODUCTION 

According to our current understanding, reionization of neutral hy- 
drogen (H I ) is most likely an extended process operating over red- 
shifts 6 < 2 < 15 and is most likely driven by stellar sources form- 
ing within galaxies (for reviews, see Choudhury & Ferrara 2006; 
Choudhury 2009). The H I distribution at these epochs is character- 
ized by ionized bubbles centered around galaxies, with character- 
istic sizes determined by the efficiency of photon production and 
the clustering properties of the galaxies (Furlanetto, Zaldarriaga & 
Hemquist, 2004). These bubbles overlap as reionization proceeds, 
and the H II (ionized hydrogen) distribution develops a complex 
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topology. Statistical probes of the redshifted 21-cm signal, like the 
power spectrum (see Furlanetto, Oh & Briggs 2006 for a review) 
and the bispectrum (Bharadwaj & Pandey, 2005), carry signatures 
of these bubbles. A detection of the statistical signal would indi- 
rectly constrain the bubble distribution. In contrast, it would be pos- 
sible to directly probe this by detecting the individual H II regions. 

At this moment it is not clear whether the first generation of 
21-cm observations would be able to detect the detailed topology 
of H II regions and constrain reionization histories. Most likely, the 
present day instruments (GMRT^ (Swarup et al., 1991) ) or those 
expected in the near future (e.g., MWA^, LOFAR^) might just be 
adequate for detecting the presence or absence of an ionized re- 
gion. The first target could possibly be detection of H II bubbles 

^ http://www.gmrt.ncra.tifr.res.in 

^ http://www.haystack.mit.edu/ast/arrays/mwa/ 

^ http://www.lofar.org/ 
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around targeted luminous sources (say, quasars) which can have 
much larger sizes than the typical galaxy generated ionized regions. 
The growth of such a region may, in the simplest situation, be mod- 
eled as a spherical ionized bubble embedded in a uniform neutral 
medium and its growth can be followed analytically (Shapiro & 
Giroux, 1987). The redshifted H 1 21-cm signal from such a bubble 
will be buried in foregrounds and noise, both of which are consid- 
erably larger than the signal. Further the noise in different pixels of 
a radio interferometric image is correlated. It is a big challenge to 
detect this faint redshifted 21-cm signal of an ionized bubble. Datta, 
Bharadwaj & Choudhury (2007) (hitherto Paper I) have developed 
a visibility based matched filter technique that optimally combines 
the H I signal of an extended bubble while removing foregrounds 
and minimizing the noise. The analysis shows that density fluctua- 
tions in the neutral hydrogen outside the bubble impose a limit on 
the comoving radius of the smallest bubble that can be detected. 
This limit is independent of the observing time. Simulations show 
(Datta, Majumdar, Bharadwaj & Choudhury 2008; hitherto Paper 
II) that bubbles of comoving radius < 6 Mpc and < 12 Mpc can- 
not respectively be detected using the GMRT and the MWA respec- 
tively, however large be the observing time. In addition, it is also 
possible to show that the redshift 2 ~ 8 is optimum for bubble de- 
tection (Datta, Bharadwaj & Choudhury 2009; hitherto Paper III). 
At the redshift 2 = 8, for a bubble located at the center of the field 
of view (FoV), it will be possible to detect (3a) a bubble of comov- 
ing radius > 19 Mpc and > 22 Mpc with 1000 hrs of observation 
using the GMRT and MWA respectively provided the gas outside 
the bubble is completely neutral. This prediction is somewhat mod- 
ified if a significant fraction of the neutral gas outside the bubble is 
ionized by 2 8. In a situation where the neutral hydrogen fraction 
outside the bubble is xh i ~ 0.6 at 2 = 8 (Choudhury & Ferrara, 
2006; Choudhury, 2009), the comoving radius of the smallest bub- 
ble that can be detected in 1000 hrs of observation is ~ 24 Mpc 
and ~ 28 Mpc for the GMRT and MWA respectively. 

Most of the analysis of matched filter bubble detection (Pa- 
pers I, II and III) assumes the ionized region to be a spherical bub- 
ble. However, a growing spherical bubble will appear anisotropic 
for a present day observer due to the finite light travel time 
(FLTT) (Wyithe & Loeb, 2004; Yu, 2005; Wyithe, Loeb & Barnes, 
2005; Shapiro et al., 2006; Sethi & Haiman, 2008), and due to 
the evolution of the global ionized fraction (Geil et al., 2008). 
These anisotropies, if detected, would provide important informa- 
tion about the quasar luminosity and the evolution of the global ion- 
ization fraction. While this is an interesting possibility, it is unlikely 
that it will be feasible to discern such details using either the present 
day instruments (GMRT) or those expected in the near future. (eg. 
MWA) which are just adequate for detecting the presence or ab- 
sence of an ionized bubble. Given this, it would be most appropri- 
ate to search for ionized regions using a spherical bubble as a tem- 
plate in the matched filter technique. The search can either be com- 
pletely blind where we vary the four parameters [Ox, dy,zt, Rb], or 
it can be targeted along a known quasar where it is sufficient to 
vary [zb , Rb] ■ Here 6x,9y, Zb are respectively the two angular coor- 
dinates and redshift of the bubble center, while Rb is the bubble's 
comoving radius. In the matched filter analysis, the signal to noise 
ratio is maximum when the parameters of the filter exactly match 
those of the bubble actually present in the data. In other words, a de- 
tection is achieved by varying the parameters of the filter so that it is 
exactly matched to the signal of the bubble that is actually present 
in the data. The apparent anisotropies will, however, introduce a 
mismatch between the signal and the filters. This is potentially a 
serious issue for both a blind search and targeted bubble detection. 
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Figure 1. The dashed lines show different stages of the growing H II bubble 
ai'ound a quasar Q. The solid curve shows the apparent shape of the 
H II bubble for a distant observer O located in the direction indicated by 
the arrows. 

In this paper we study the impact of this apparent anisotropy on 
our ability to detect ionized bubbles in a targeted search around a 
known quasar using the matched filter technique. 

In reality apart from these anisotropies the ionized regions 
around quasars would not be spherical because of various other 
effects like inhomogeneities in the IGM, overlapping ionized bub- 
bles arising from stellar sources and possible anisotropic emission 
from the quasar itself. Moreover, the foreground subtraction could 
also subtract a part of the signal, which may change the recovered 
bubble shape. The anisotropy in bubble shape caused by these ef- 
fects would vary from bubble to bubble, and it is rather difficult to 
incorporate these effects in our analytical estimations, so we do not 
attempt this here. Geil & Wyithe (2008) and Geil et al. (2008) have 
studied several of these effects using semi-numerical simulations. 

A brief outline of the paper follows. In Section 2. we discuss 
the equations governing the growth of a spherical ionized bubble 
and the apparent anisotropy due to the FLTT. Section 3. reviews 
the matched filter technique for bubble detection. In Section 4. we 
quantify the apparent anisotropy due to the FLTT and assess its im- 
pact on bubble detection. We also discuss the effect of an evolving 
global neutral fraction. We summarize and conclude in section 5. 
Growth equation for ionized bubble with an evolving neutral frac- 
tion has been discussed in Appendix A. 

Unless mentioned otherwise, throughout the paper we present 
results for the redshift 2 = 8 assuming the neutral hydrogen frac- 
tion Xh I ~ 0.6 outside the bubble, with the cosmological param- 
eters h = 0.74, Qm = 0.3, = 0.7, Qbh'^ = 0.0223. 



2 GROWTH OF IONIZED BUBBLES AND THEIR 
APPARENT SHAPE 

We first review the equation governing the growth of a spheri- 
cal H II region around a quasar (Shapiro & Giroux, 1987; White, 
Becker, Fan & Strauss, 2003; Wyithe & Loeb, 2004; Wyithe, Loeb 
& Barnes, 2005; Yu, 2005; Shapiro et al., 2006; Sethi & Haiman, 
2008). For our purpose we consider the form given in eq. (7) of Yu 
2005, which is 

4:7V d / 1 \3\ f \ ^ ^/ \23 

"3"^ (a^Hi [nH)r ) = iVphslr) - -naBC[nH) r . (1) 

Here the quasar is assumed to be triggered at a cosmic time ti, 
and T = t — ti denotes the quasar's age at any later time t. The 
variable r(r) denotes the radius of the spherical ionizing front 
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Figure 2. The growth of an H II bubble for Models (i) and (ii) shown by the 
solid and dashed curves respectively. 



(Figure 1) at the instant when a photon that was emitted from the 
quasar at r catches up with the ionizing front. The quasar is as- 
sumed to emit ionizing photons isotropically at a rate Nphs{T) 
which completely ionizes the hydrogen inside the spherical bub- 
ble of radius r'(r). The bubble is surrounded by gas with mean hy- 
drogen number density {uh) and neutral fraction xh i • The term 
as(= 2.6 X 10^^'^ cm^ s~^) is the recombination coefficient to 
excited levels of hydrogen at T = 10"* K, and C = {rijj , ) / {nH)"^ 
is the clumping factor, which quantifies the effective dumpiness 
of the hydrogen inside the bubble. Eq. (1) essentially tells us that 
the growth of the bubble is driven by the supply of ionizing pho- 
tons after accounting for the photons required to compensate for 
the recombinations inside the existing ionized region. The effect of 
Hubble expansion is not included in this equation. 

Assuming that all the quantities except r(r) and Nphsir) are 
constant over the time-scale of the bubble's growth, we have the 
solution 



Nphsir) exp 



4:Tr{nH)xH I 



where T,.ec is the recombination time defined as 

Tree = Xhi [C {uh) OLb)^^ 

^ -.xio-,.(25L) 



dr 



(2) 



7.4 

l + z 



(3) 



For the redshift of our interest {z — 8) and assuming that xh i = 
0.6 and C = 30 (Yu & Lu, 2005) we have Tree ^ 1.33 x lO'^yr 
which we use throughout. 

Following Yu 2005, we consider two different models for 
the photon emission rate NphaiT). In Model (i) NphaiT) = 
Nphs,i is a constant whereas it increases exponentially in Model 
(ii) Nphsir) = Nphs,i exp{T/Ts). Here rs ~ 4.5 x 
10^ yr [^ Q i(i-g) ] whsre rg is the Salpeter time-Scale and e is 
the mass to energy conversion efficiency whose value typically is 
e ~ 0.1 which we adopt throughout (Yu & Tremaine, 2002; Yu & 
Lu, 2004, 2005). This gives Tb — 5 x 10^ yr which we have used 
throughout in Model {ii). The bubble growth r(T) for model (i) 
and (ii) are respectively given by 



r(r) = rs 



exp 



(4) 
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Figure 3. The variation of r with 
(ii) (right)) with Ts = 88 Mpc. In both panels tq takes values 
(1.3,2.7,5.3,8.0,10.6, 13.3) x 10^ yr from the bottom to top. The values 
of Ts and r are both in comoving Mpc. 
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Figure 4. The shape of the H II bubble as it would appear in a redshifted 
21 -cm map. The distant observer is located to the right of the figure, and 
the line of sight is along the x axis. The quasar is located at the intersection 
of the two dashed lines. The photon emission model, rs and tq values are 
the same as in Figure 3. The left and right panels correspond to Model (i) 
and (ii) respectively. 



and 

r(r) = rs 
where 



exp 1 — ) — exp ( — 

'^s / \ 'Tre.c 



rs = 



3Nphs ,i 'Tr, 



4:1T X} 



{uh) 



(5) 



(6) 



for which the solutions are shown in Figure 2. In both mod- 
els the bubble has an initial period (r < Tree) of rapid growth 
r(T) ~ [t /TrecY^'^ ■ In Modcl (i) the bubble radius subsequently 
approaches a constant value r(r) = r^. In Model [ii) the bub- 
ble radius does not reach a steady value but continues to grow at a 
slower rate beyond r > . 

To visualize the apparent shape of the bubble as seen by a 
present day observer we need the relation between r and <j), where 
<j) is the angle between observer's line of sight (LOS) and the point 
A under consideration on the ionization front (Figure 1). The light 
travel time starting from the quasar at r to the point A and then to 
the present day observer is [r (r ) /c] ( 1 — cos (j>) more compared to 
the photon that was emitted from the quasar at age tq and travels 
straight to the present day observer. This gives 



TQ = T + 



(7) 



(eq. (3) of Yu 2005) which we use along with eq. (4) or eq. (5) to 
determine r as a function of (f> for Models (i) and (ii) respectively. 

The photons from the front part of the bubble {(f) < 90°) take 
less time to reach the observer as compared to photons from the rear 
part (0 > 90°). As a consequence the observer sees different parts 
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Figure 5. This shows the apparent shape of an H 11 bubble as it would to a 
distant observer located along the line of sight (LOS) indicated in the figure. 
The apparent center of the bubble shifts by QC = Ar/2 from the quasar 
Q to the point C. The apparent radius is CB = along the LOS and 
CA = R± in the perpendicular direction. 



of the bubble at different stages of its growth whereby it appears 
anisotropic along the LOS. 

The apparent shape of the H II region is controlled by the 
parameters and tq, where tq is the quasar's age as seen by 
a present day observer. In Figure 3 we have explicitly shown r 
as a function of </> for rs — 8 x crrec(l + z) — 88Mpc 
(comoving) and Tg/rrec = 1, 2, 4, 6, 8, 10, which correspond to 
TQ = (1.3, 2.7, 5.3, 8.0, 10.6, 13.3) x lO'^ yr respectively In Fig- 
ure 4 we show a radial section through the center of the bubble as 
it would appear in a redshifted 21-cm map. The difference in light 
travel time across the bubble is more for larger bubbles, which is 
why we have chosen a particularly large value of rs to illustrate the 
anisotropy. Note that here and in all subsequent discussion we use 
comoving length-scales. 

The bubble grows rapidly when Tg/Trec < 1 which is close to 
the instant when the quasar was triggered, We then have a large dif- 
ference between the front and back surfaces (Figure 3). The back 
surface, is viewed at an earlier phase of growth compared to the 
front surface. This has mainly two effects on the bubble's apparent 
shape (Figure 4), (i) the center is shifted along the LOS towards 
the observer, (ii) the bubble appears anisotropic. We also note that 
the bubble appears elongated along the LOS in the early stage of 
growth. The difference between the front and back surfaces gradu- 
ally comes down with increasing tq. The apparent center then ap- 
proaches the quasar position, and the elongation along the LOS also 
diminishes. While Models (i) and (ii) are nearly indistinguishable 
at small tq, the behavior are somewhat different at large tq. 

For a fixed rs and tq, we use rmin and rmax to denote the 
bubble's comoving radius r at = 180° and 0° respectively. The 
difference Ar — rmax — rmin causes the bubble's center to shift 
by a comoving distance Ar/2 relative to the quasar along the LOS 
towards the observer (Figure 5). The bubble has comoving radii 
R\\ = {vmax + rmin)/^ and R± along the LOS and perpendicular 
to it respectively (Figure 5). In the subsequent discussion we use 
R± to characterize the comoving size of the anisotropic bubble. 
The dimensionless ratio 
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Figure 6. Contours of R± (dashed lines) and R f (soUd lines) in comoving 
Mpc, shown for Models (i) and (ii) in the bottom and top panels respec- 
tively. R.f (defined in Section 4) is the radius of the spherical bubble that 
best matches an anisotropic bubble with parameters rs and tq . 



R± . Further, we use 



R± 



(9) 



2R± 



expresses the shift in the bubble's center as a fraction of its radius 



to quantify the anisotropy of the bubble. A value > indicates 
that the bubble appears elongated along the LOS, whereas < 
indicates that it appears compressed along the LOS. We have calcu- 
lated R±,s and rj for a range of r s and tq , for which the results are 
shown in the contour plots of Figures 6, 7 and 8 respectively. The 
rs, TQ range was chosen so that the largest bubble has a comoving 
radius roughly in the range 70 — 80 Mpc. 

In our models the bubble, as viewed in the quasar's rest-frame, 
has a short period of rapid growth r < 10^ yr beyond which it sat- 
urates (Model (i)) or grows slowly (Model (ii)) (Figure 2). Naively 
one would expect the bubble's image in redshifted 21-cm maps to 
exhibit a similar behavior. However, we see that the bubble's image 
continues to grow (Figure 6) well beyond r = 10^ yr. The effect 
is more for larger bubbles where the light travel time is longer. The 
shift and anisotropy also are non-zero well beyond the period of the 
bubble's actual growth. For a fixed r^, the shift decreases monoton- 
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Figure 7. The dashed Unes show contours of the shift parameter s defined 
by Eq. (8). The solid lines show contours of So which quantifies the shift 
between the center of the best matched filter and the quasar (defined in 
Section 4.). Models (i) and (ij) are shown in the bottom and top panels 
respectively. 



ically with increasing tq. While the shift can be as large as s > 1 
(Figure 7) for a large bubble seen in an early phase of its growth, 
the typically value of s is in the range 0.2 < s < 0.6. The bubbles 
appear elongated ( > 0). in the early phase of their growth (Fig- 
ure 8). For a fixed tq, we have more anisotropy for a larger bubble. 
The elongation diminishes with increasing tq . This is followed by 
a situation where the bubble appears compressed (77 < 0) along the 
LOS. This occurs when a bubble is viewed in the late stage of its 
growth where most of the bubble, except a small part of the back 
surface, has nearly stopped growing. The transition from an elon- 
gated bubble to a compressed bubbles is clearly visible in Figure 4. 
For both Models (i) and (ii), by and large, we expect the bubbles 
to be compressed with anisotropies of the order of | r\ |~ 0.1 in 
Model (i), and somewhat larger anisotropies in Model (ii). Elon- 
gated bubbles may have anisotropies of the order ~ 0.1 — 0.5, 
but these will be seen only in bubbles surrounding recently trig- 
gered quasars. 
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Figure 8. This shows contours of the anisotropy parameter r\ (defined by 
Eq. (9)) for Models (i) and (ii) in the bottom and top panels respectively. 
The contours are at an interval of 0.1. 



3 MATCHED FILTER BUBBLE DETECTION 
TECHNIQUE 

The quantity measured in radio-interferometric observations is the 
visibility V((7, v) which is related to the specific intensity pattern 
on the sky Iv(&) as 



V{U,v) = 



(10) 



Here the baseline U = d/A denotes the antenna separation d 
projected in the plane perpendicular to the line of sight in units 
of the observing wavelength A, S is a two dimensional vector in 
the plane of the sky with origin at the center of the FoV (i.e. 
phase center), and Aid) is the beam pattern of the individual an- 
tenna. For the GMRT this can be well approximated by Gaussian 
A% = 6-*'/"°' where ^ 0.6 eipwHTvi- 

The visibility recorded in a radio-interferometric observations 
is actually a combination of several contributions 

V(Jj,v) = S{U,v) + HF(U,v) + N{U,v) + F{U,i^). (11) 

where the first term S(U, u) is the expected signal of the ionized 
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region that we are trying to detect. Note tliat tlie signal from an 
ionized region appears as a decrement witli respect to the uniform 
background 21 -cm radiation. This uniform background is however 
recorded only by the baseline at zero spacing which is usually not 
considered in radio-interferometric observations. Readers are re- 
ferred to Paper I (section 2) for more detail regarding how to calcu- 
late S{U, v). HF{U, v) is the contribution from fluctuations in the 
H I distribution outside the ionized bubble, N{U , u) and F{U , u) 
are the noise and foreground contributions respectively. 

The signal from an ionized bubble is expected to be buried 
deep under the other contributions (noise and foreground) which 
typically are orders of magnitude larger. We introduce an estima- 
tor E[6x, 6y,Zfy, Rfy] to scarch if a particular ionized bubble (with 
bubble parameters [9x, Oy,Zb, Rb\ ) is present in our observation. 
The estimator combines all the visibilities weighted by a filter 



(12) 



Note that in order to keep the notation compact we do not explic- 
itly show the parameters [9x, 6y,zi,, Rt] and V{Ua,i^b) represents 
the observed visibility. The filter is designed (Paper I) to optimally 
combine the signal corresponding the the bubble that we are trying 
to detect while minimizing the contribution from the other contam- 
inants. The expectation value of E is 



(E) =Y.^HUa,lyb)S{tja,lyb) 



(13) 



The other contributions to V are assumed to be random variables 
of zero mean, uncorrelated to the filter, and hence they contribute 
only to the variance 

((A£;)') = {{AEf)HF + ((A£)')iv + {iAEf)F (14) 

where the subscripts HF, N, F respectively refer to the contribu- 
tions from the H I fluctuations, noise and foregrounds. The signal 
to noise ratio for the estimator is defined as 



SNR: 



(E) 



{(Ai5)2> 



(15) 



Bubble detection will be carried out by analyzing the SNR for 
filters with different values of the parameters [9x,9y, zt, Rb] - If an 
ionized bubble is actually present in the FoV, the SNR will peak 
when the filter parameters exactly match the parameters of the bub- 
ble. We shall have a statistically significant bubble detection if the 
peak SNR > 3 or SNR > 5 for a 3cr or 5a detection respectively. 

While the filter assumes the bubble to be spherical, it is ev- 
ident from the previous section that bubbles are expected to ap- 
pear anisotropic in redshifted 21-cm maps. It is in principle possi- 
ble to introduce the anisotropy as an additional search parameter. 
This, however, would be required for preliminary bubble detection 
only if there is a severe mismatch between the spherical filter and 
the anisotropic image resulting in considerable degradation of the 
SNR. 



4 RESULTS 

In order to keep the computational time of the present investiga- 
tion within bounds we consider a targeted search along the LOS 
to a known high z quasar The FLTT does not affect the angular 
position of the bubble's center and 9x,9y will coincide with the 
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Figure 9. This shows anisotropic shapes and corresponding expected SNR 
contours for Model (i) considering 1000 hrs of observation with GMRT. 
From bottom to top solid curves in the left panels correspond to different 
apparent shapes of ionized regions for tq = (1.3, 5.3, 10.6) X 10*^ yr re- 
spectively with Ts = 88 Mpc. The dotted curves represent spherical filters 
for which the SNR is maximum for the shape in that con'esponding panel. 
Intersecting point of the dash-dotted lines represent position of the quasar 
and intersecting point of the dash-dotted and solid lines represent position 
of the filter center Right panels show the matched filter SNR contours for 
corresponding left panel anisotropic shapes as function of filter radius R f 
and filter center frequency coordinate v j . 



angular position of the quasar. We consider a search varying the 
parameters Rf and i/f which respectively correspond to the filter 
radius and the position of the filter's center along the LOS. In the 
bottom panels of Figure 9 we have explicitly shown the the result 
for a situation where a bubble with parameters rs = 88 Mpc and 
Tq — 1.3 X lO'^yr is present surrounding a quasar which is lo- 
cated at vq — 158 MHz in the center of the FoV. The right panel 
shows the SNR for a search in the parameters and Rf. We find 
that the SNR peaks at lyj = 160.4 MHz and Rf = 32 Mpc, with 
SNR = 3.99 for 1000 hr of observation. The shape of the best 
match filter is shown in the left panel along with the bubble's appar- 
ent shape. The middle and top panels consider the same bubble at 
two later stages of its growth. We see that the bubble, which appears 
elongated (rj = 0.48) in the bottom panel, becomes compressed in 
the middle and top panels (rj = —0.11, —0.2 respectively). In all 
cases the best matched filter provides a reasonably good represen- 
tation of the bubble that is actually present. The mismatch between 
the filter and the bubble that is actually present causes the SNR to 
be lower than that which is expected if the image of the bubble were 
spherical. 

We have considered bubble detection for a range of rs and tq 
values for which the radius of the best match filter Rf are shown 
in Figure 6. Figure 7 shows the observed shift So = (rq — rc)/Rf 
where rg and are the comoving distances to the quasar and the 
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Figure 10. The dashed contours show / (defined in Section 4.), while the 
solid contours demarcate the regions where a Scr and 5<j detection is possi- 
ble. The bottom and top panels show results for Models («) and {ii) respec- 
tively. 



center of the filter respectively. We see that for most of the parame- 
ter space Rf and So for the best match filter closely follow R± and 
s of the corresponding bubble image. 

The situation where both the bubble and the corresponding 
quasar have been detected presents an interesting possibility. Both 
Rf and So are known in such a situation. The curves corresponding 
to constant Rf and So are non-degenerate, and it is possible to lo- 
cate an unique point in the Ts — tq parameter space once Rf and So 
are known. This holds the possibility of allowing us to estimate the 
age of the quasar responsible for the bubble. While the age estimate 
would depend on the model assumed for the photon emission rate, 
our analysis considering two different models indicates that we do 
not expect this variation to be very large. 

The mismatch between the anisotropic bubble and the spher- 
ical filter is expected to cause a reduction in the SNR. For the 
best match filter Figure 10 shows the ratio / of the SNR for an 
anisotropic bubble to the SNR for a spherical bubble with the same 
Rf. We find that this ratio has value in the range / ~ 0.9 — 1.0 for 
nearly the entire parameter space that we have considered. There 
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Figure 11. This shows the evolution of xu i across the redshift range that 
is relevant for the growth and detection of a bubble with zq = 8. The lower 
and upper dotted curves respectively show the early and late reionization 
models proposed by Choudhury & Fen'ara (2005, 2006). The solid curve 
shows the model adopted for our analysis. 

is a small region with small tq, which corresponds to rapid bubble 
growth, where the ratio is in the range / ~ 0.8 — 0.9. We note that 
this matches with the region where the anisotropy is large i] > 0.4. 
For 1000 hr of observations with the GMRT, the Vs, tq range for 
3(7 and 5cr detection are shown in Figure 10. The smallest filter 
radius (comoving) for a 3cr and a 5cr detection are approximately 
24 Mpc and 33 Mpc respectively, which matches with the corre- 
sponding values when the finite light travel time is not taken into 
account. 



4.1 Anisotropy due to evolving xh i 

The entire analysis, till now, assumes a fixed neutral fraction 
Xh I ~ 0.6 outside the bubble. The apparent anisotropy of a spher- 
ical bubble will, in principle, be affected by the evolution of xh i ■ 
The bubble's front surface which is seen at a lower redshift will 
have a lower neutral fraction compared to the back surface. This 
will further increase the difference between the radius of the front 
and back surfaces as compared to the situation where xh i ~ 0.6 
for which the results have been shown in Figure 3 and 4. 

The evolution of xh i during reionization is largely unknown, 
and there exists a wide range of evolution histories all of which 
are consistent with the presently available observational constrains. 
The lower and upper dotted curves in Figure 1 1 respectively shows 
the early and late reionization models proposed by Choudhury 
& Ferrara (2005, 2006). The entire region bounded by these two 
curves is consistent with the current observational constraints. To 
assess the impact of xh i evolution on bubble detection we have 
considered a simple model (the solid curve in Figure 11) where 
xh I ^ evolves linearly in the vicinity of z = 8. We do not expect 
Xh I to decline very rapidly around 2 = 8, and the slope has been 
chosen so as to approximately match those of the two dotted curves 
in Figure 1 1 . 

The subsequent discussion of this subsection is entirely re- 
stricted to the largest bubble (rs — 100 Mpc) for which the bub- 
ble's comoving radius grows to 100 Mpc by r = 10* yr when 
Xh I is constant at 0.6. The bubble's apparent image spans from 
z = 8.2 to 7.6, and xh i changes from 0.64 to 0.52 across the 
bubble in the evolving xh i model that we have considered. For 
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Figure 12. This dotted curve shows the growth of an H II bubble for our 
model of evolving global neutral fraction. The quasar age is set to tq = 
1.3 X 10^ yr at 2 = 8. The results for a constant neutral fraction xh i = 
0.6 are shown for compaiison (solid curve). The bottom and top panels 
respectively correspond to Models (i) and (it). 
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Figure 13. This shows the variation of r] (top panel) and s (bottom panel) 
with TQ for an H II bubble with Vs = 100 Mpc. In both panels the solid 
and dashed curves represent the constant and evolving xh i models respec- 
tively. The two different sets of curves in each panel coiTespond to Models 
(i) and (ii) respectively. 



the same model, xh i is 0.7 at the redshift where the bubble is 
born if the quasar has an age of tq — 1.3 x 10* yr at 2 = 8. 

We have modified eq. (1) to incorporate the evolution of xh i 
(see Appendix A for details) for which Figure 12 shows the growth 
of the bubble's comoving radius. For both Models (i) and (ii) we 
find that incorporating the evolution of xh i initially slows the 
growth of the bubble relative to the situation where xh i is con- 
stant. This is reversed in the later stages where the growth is en- 
hanced with respect to the situation where xh i is constant. For 
the parameter range that we have considered, the difference in the 
bubble radius due to the evolution of xh i is found to be less than 
5%. Further, there is a significant difference only at the later stages 
of the bubble's growth (r > 3r,,ec ~ 4 x 10^ yr). We next con- 
sider the effect of the evolution of xh i on the anisotropy and shift 
parameters ri and s (Figure 13). We find that the effect is maximum 
when the bubble is seen in the late stage of its growth where the ap- 
parent image is compressed along the line of sight. The evolution 
of Xh I acts to reduce the compression and increase the shift rela- 



tive to a situation where xh i is constant. The maximum change is 
found to be around 10%. This change in the anisotropy has a very 
small impact on matched filter bubble detection. The radius of the 
best matched filter, we find, changes by less than 5% whereas the 
corresponding SNR changes by less than 1%. The changes being 
extremely small, we have not explicitly shown these results. 

The discussion till now has focused on a large bubble with 
Ts = 100 Mpc. The effect of the evolution of xh i will be less for 
smaller bubbles. 



5 SUMMARY & CONCLUSIONS 

The apparent image of a quasar generated, spherical ionized re- 
gion appears anisotropic due to the finite light travel time and the 
evolution of the global neutral fraction. We have quantified this 
anisotropy and assessed its impact on matched filter bubble detec- 
tion in redshifted 21 -cm maps at 2 = 8. We find that the anisotropy 
is less than 20% for most of the parameter space that we have con- 
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sidered. Further, the bubble appears compressed along the LOS for 
nearly the entire parameter range where a detection is possible. In 
addition to the anisotropy, the bubble's center is also shifted from 
the actual position of the quasar along the LOS. This shift, which 
can be rather large and may exceed the bubble radius, typically 
varies across the range 20 — 80%. We find that these effects do not 
cause a significant degradation of the SNR in matched filter bubble 
search and spherical filters are adequate for a preliminary bubble 
detection. We show that such a detection where the radius of the 
best match filter, its center and the quasar redshift are all known 
can be used to estimate tq the age of the quasar. This detection 
also provides an estimate of the neutral fraction xh i provided the 
quasar luminosity is known. Follow up observations, subsequent 
to a detection, can in principle be used to measure the anisotropy 
which would impose further constrains on the quasar luminosity, 
TQ and Xh i . The global neutral fraction is not expected to evolve 
very rapidly with z around 2 = 8, and this make a rather small 
change (< 10%) for the anisotropy and the shift relative a situa- 
tion where xh i is constant. 

An earlier work (Paper III) has shown that both the GMRT or 
the MWA are most sensitive to the signal from an ionized bubble 
if it is at a redshift z 8. At this z a bubble of comoving radius 

> 24 Mpc and > 33 Mpc can respectively be detected at Scr and 
5(T significance with 1000 hrs of observation. This, however, pre- 
supposes the presence of an ionized bubble located near the center 
of the FoV. It is, however, necessary to take into account the abun- 
dance of ionized bubbles when considering possible observational 
strategies. 

How many H II bubbles do we expect to find above the detec- 
tion threshold in a single GMRT observational volume at redshift 
2 = 8? There is, till date, no quasar detected with z > 6.4, though 
there are galaxies which have been detected at 2 > 8 (eg. Lehnert et 
al. (2010)) and it is quite likely that quasars may be detected in the 
near future. At present one thus typically uses theoretical models 
which correctly predict the observed quasar abundance and evolu- 
tion at 2 < 6 to predict the expected number of detectable quasar 
H II bubbles at higher redshifts, though these predictions are highly 
model dependent and uncertain (Rhook & Haehnelt, 2006). Pre- 
dictions based on the CDM merger driven model (Wyithe & Loeb, 
2003; Wyithe, Loeb & Barnes, 2005; Rhook & Haehnelt, 2006) 
indicate that we expect 1 ionized bubble of comoving radius 

> 20 Mpc around an active quasar in 400 degree^ and 16 MHz 
volume. This means that the probability of a detection is less than 
5% in a single random GMRT observation volume (~ 10 degree^ 
and 32 MHz) at redshift 2 = 8. Including the possibility of fossil 
bubbles whose quasars have switched off increases the number of 
H II bubbles (Wyithe, Loeb & Barnes, 2005), and we expect ^ 5 
detectable bubbles in a random GMRT observation. Note that, the 
shape of these fossil bubbles would not be distorted due to the FLTT 
as they have ceased to grow. 

The typical size of quasar generated H II bubbles at 2 ~ 8 is 
an important factor in matched filter bubble detection. Maselli et al. 
(2007) predict the comoving radius of quasar generated H II regions 
to be ~ 45 Mpc at 2 = 6.1 with xh i =0.1. This size is expected 
to be less at 2 ~ 8 where the global neutral fraction is higher. The 
typical radius of a galaxy generated bubble is expected to be ~ 10 
comoving Mpc at 2 ~ 8 (McQuinn et al., 2007). These galaxy gen- 
erated bubbles outside the quasar generated bubble that is targeted 
for detection will contribute to the fluctuations in the matched fil- 
ter estimator E. However, simulations (Paper II) indicate that small 
bubbles of size 10 Mpc will not have much impact on detection in 
the range (> 24 Mpc) where a detection is feasible. 



Finally, we note that there are various other sources of 
anisotropy such as the spatial fluctuations in the neutral hydro- 
gen distribution due to galaxy generated generated ionized bubbles, 
anisotropic emission from the quasar itself which have been ig- 
nored, galaxy generated generated ionized bubbles in the surround- 
ing IGM can cause significant anisotropies in the quasar generated 
ionized bubbles. One needs to use detailed simulation to investigate 
this which is beyond the scope of this paper. We plan to address 
these issues in future. 
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APPENDIX A: GROWTH EQUATION FOR H ii BUBBLE 
CONSIDERING EVOLVING NEUTRAL FRACTION 

Here we consider the growth of an H II bubble embedded in a 
medium whose neutral fraction xh i outside the bubble is evolv- 
ing. The evolution of the neutral fraction is parameterized using 



XHl (r) = XHl {tq) f(T) 



(Al) 



where xh i (tq) is the neutral fraction at the redshift zq — 8, and 
the function /(r) which quantifies the evolution satisfies fijo) = 
1. Further, the ionizing photon emission rate of the quasar is ex- 
pressed as Nphs ~ Nphs,i s(r) with s(r) — 1 for Model (i) and 
s(r) = exp{T/Ts) for Model (ii). Using these, the equation gov- 
erning the growth of an ionized bubble (eq. 1) may be written as 



d_ 

d7 



(A2) 



where Tree and are constants calculated using eq. (3) and eq. (6) 
respectively using xh i (tq) i.e. value of xh i at 2 = 8. Note that 
we have ignored the effect of the finite light travel time from the 
quasar to the ionization front on /(r). We then have the solution 



exp 



1 

Tree 



dr" 
fir") 



sir') 



^^'^^ ^" ^ ./(t) / Tree 


(A3) 

The analysis is considerably simplified if we approximate the evo- 
lution of the neutral fraction as 



1 + b 



T - TQ 



fir) \ Tree 

whereby eq. (A3) can be written in terms of a single integral as 



(A4) 
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1 - b 



TQ 



b / T — t' 

2 V Tree 



s(r)^ (A5) 

Tree I 



The integral in eq. (A5) can be easily evaluated numerically. We 
use this equation to study the growth of an H II bubble for a wide 
range of values of the parameters Vs and tq , taking into account the 
variation of the neutral fraction. It is quite evident that we recover 
eqs. (4) and (5) in the situation where 6 = i.e. the neutral fraction 
does not evolve. 

The evolution of the neutral fraction xh i is largely unknown 
at the redshifts 2 ~ 8. For the purpose of our discussion we con- 
sider two models proposed by Choudhury & Femra (2005, 2006) 
which implements most of the relevant physics governing the ther- 
mal and ionization history of the IGM through a semi-analytical 
formalism. Figure 11 shows the predicted evolution of xh i over 
the relevant z range. For the subsequent analysis in this paper we 
have used 6 = 0.015 in eq. (A4) whereby the predicted xh i (2) is 
in rough consistency with the evolution allowed by the two models. 

The growth of an H II bubble with evolving xh i (dotted line) 
has been shown in Figure 12 for both Model (i) and (ii) of photon 
emission rate. For comparison with our previous analysis in this 
paper, growth for constant xh i (solid line) is also shown in the 
same figure. 



